Sparse approximate matrix-matrix multiplication for density matrix purification with error control

We propose an accelerated density matrix purification scheme with error control. The method makes use of the scale-and-fold acceleration technique and screening submatrix products in block-sparse matrix-matrix multiplies to reduce computational cost. An bound a parameter sweep are combined select threshold value for screening, such that can be controlled. evaluate performance comparison without product screening.